Origin of the Superconductivity in /3'-(BEDT-TTF) 2 ICl2 under High Pressure 
and /3-(BEDT-TTF) 2 X at Ambient Pressure 
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We present a theoretical study on the superconductivity in /3'-(BEDT-TTF) 2 ICl 2 at T c = 14.2K 
under high hydrostatic pressure found by Taniguchi et al., which is the highest record among organic 
superconductors. Its electronic structure is well expressed by the anisotropic triangular lattice 
Hubbard model at half filling. In the present work, we study this effective model by using the 
fluctuation-exchange (FLEX) approximation. In the obtained phase diagram, the superconductivity 
with d x 2- y 2 like symmetry is realized next to the antiferromagnetic (AF) insulating phase, as a result 
of the 1D-2D dimensional crossover driven by the pressure. The obtained maximum T c is 16~18K. 
In addition, the superconductivity in /3-(BEDT-TTF) 2 X is also understood in the same framework. 

PACS numbers: 74.70.Kn, 74.70.-b 
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Up to now, many kinds of organic superconductors 
composed of BEDT-TTF (abbreviated as ET) molecules 
have been discovered and studied intensively. Especially, 
k-(ET) 2 X compounds attract much attention because 
it shows the maximum superconducting (SC) transition 
temperature (T c rj 12K) among them [1]. Recently, the 
temperature-pressure phase diagram of k-(ET) 2 X salt 
could be produced by the fluctuation-exchange (FLEX) 
approximation [2-4]. It is a kind of the self-consistent 
spin-fluctuation theory [5,6]. According to the FLEX ap- 
proximation, d-wave superconductivity is expected which 
is mediated by the strong antiferromagnetic (AF) fluctu- 
ations due to the Coulomb interaction. 

Quite recently, Taniguchi et al. found the supercon- 
ductivity at T c = 14.2K in /3'-(ET) 2 ICl 2 under high hy- 
drostatic pressure (P > 8.2GPa), which established the 
new record of T c among the organic superconductors [7] . 
At ambient pressure, /3'-(ET) 2 ICl 2 shows a semiconduct- 
ing conductivity below the room temperature, and indi- 
cates a magnetic transition at Tn=22K. As the pressure 
is applied, the resistivity decreases gradually, and the 
metallic behavior (dp/dT > 0) is observed above Tmit 
under 6.5GPa. Note that Tmit, which is a crossover tem- 
perature, is expected to be higher than Tn- At 8.2GPa, 
insulating phase (or Tmit) disappears and the SC tran- 
sition occurs at T c = 14. 2K at the same time. 

According to the band calculation, the Fermi surface 
(FS) of /3'-(ET) 2 X (X=IC1 2 , BrICI, AuCl 2 ) compound is 
quasi one-dimensional (Q1D) [8,9]. It is contrastive that 
the FS of /3-(ET) 2 X is almost round, which shows the 
superconducting transition at T c = 2.7K for X=IBr 2 , at 
T c = 3.8K for X=AuI 2 and at T c = 1.5K (or 8K) for 
X=1.3 . The difference of the crystal structure between (3 
and /3' compounds arises from the fact that the anion size 
in /?' compounds is too small to retain the /3-type struc- 
ture [8]; in fact, the bond lengths for I-Cl and Au-Cl are 
2.52A and 2.27A respectively, which are smaller than the 
I-Br bond length in /3-(ET) 2 IBr 2 , 2.58A [8,9]. 

Until now, the deformation of the structure of /?'- 



(ET) 2 IC1 2 at P = 8.2GPa has not been determined ex- 
perimentally. In the present study, however, we natu- 
rally assume that its structure (gradually) approaches to 
that of /3-(ET) 2 X as the cell volume decreases under hy- 
drostatic high pressure, because the anion size becomes 
larger compared with the cell length. Note that the shape 
of the unit cells for f3- and that for /3'-structures are al- 
most same. 
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FIG. 1. Microscopic structures for /3'-(ET) 2 ICl 2 and for 
/3-(ET) 2 l3. They are well approximated as the dimer Hub- 
bard model (to, ti, t2, U e s) at half filling. In the present 
simplified figure, /3'-type structure is given by rotating each 
dimer in /3-type structure slightly in a clockwise direction. 

In the present work, we study the origin of the super- 
conductivity in /3'-(ET) 2 ICl 2 on the assumption that the 
crystal structure approaches to /3-type under hydrostatic 
pressure. We use the FLEX approximation by noticing 
the fact that the SC phase appears by destroying the AF 
phase by pressure. Our theory can reproduce a reason- 
able value of T c in /3'-(ET) 2 ICl 2 under high pressure, as 
well as a low T c in /3-(ET) 2 X at ambient pressure. 

The schematic structures of /?'-(ET) 2 X and ^-(ET) 2 X 
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are shown in Fig. 1. Each ellipse represents the ET 
molecule, and each highest occupied molecular orbital 
(HOMO) of the ET molecule possesses 1.5 electrons on 
average. The hopping parameters obtained by the band 
calculation are given in Table I: The hopping parameters 
for /?'-(ET) 2 AuCl 2 and those for /3'-(ET) 2 ICl 2 are simi- 
lar. Whereas, they are very different from the parameters 
for /3-(ET) 2 AuCl 2 although their crystal structures dif- 
fer slightly. This comes from the fact that the overlap 
integral between two ET molecules is very sensitive to 
the angle between them ((f)): It takes the local maximum 
(positive) values for <f> — ° ,60 ° and the local minimum 
(negative) values for <\> = 30 ° ,90 ° [10]. 

As a good approximation, we take account of only the 
anti-bonding orbit of each pair of the ET molecules con- 
nected by t(pl), which is more than 2.5 times larger than 
other hopping parameters [11]. Then, the original sys- 
tems is mapped onto the "dimcr model" at half filling, 
which is the anisotropic triangular lattice Hubbard model 
with four parameters (to, t\, t 2 ,(7eff) as shown in Fig. 1. 
Note that the x-y coordinate in the present model does 
not correspond to that for the original systems. The cor- 
responding hopping parameters for /3'-(ET) 2 ICl 2 and for 
/3-(ET) 2 I 3 are given in Table II. 





/3'-(ET) 2 ICl 2 


/3'-(ET) 2 AuCl 2 


/3-(ET) 2 I 3 


i(pl) 


-0.272 (cV) 


-0.264 (cV) 


-0.245 (cV) 


i(p2) 


0.016 


0.020 


-0.084 


t(ql) 


-0.100 


-0.100 


-0.127 


t(q2) 


-0.066 


-0.065 


-0.068 


t(c) 


-0.016 


-0.023 


0.050 



Table I: Hopping integrals for each systems. They are given 
by the well known empirical relation t = ~10S (eV), where 
S is the overlap integral given by the band calculation [8,9]. 





x = [/3'-(ET) 2 ICl 2 j 


x=\ [/3-(ET) 2 I 3 ] 


to 


0.050 (eV) 


0.064 (cV) 


tl 


0.017 


0.084 


*2 


-0.008 


0.042 



Table II: Hopping integrals in the dimer model for x 
x = 1, respectively. 
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Although the original on-site Coulomb interaction on a 
ET molecule is ~ lcV, the effective Coulomb interaction 
on a dimcr, U e g, is limited to ~ 2|t(pl)| ~ 0.5eV [11]. 

Hereafter, we study the dimer model given in Fig. 1 
at half filling. To analyze the pressure effect on the elec- 
tronic states, we put the hopping parameters of the dimcr 
model (to, t\, t 2 ) as follows: 



to = 00.5 + 0.014x (eV), 
ti = 00.17 + 0.067a; (eV), 
t 2 = -0.008 + 0.05x (eV), 



(1) 



where < x < 1 is a free parameter, x = (x = 1) 
corresponds to /3'-(ET) 2 ICl 2 (/3-(ET) 2 I 3 ); see table II. 



Non-linear terms with respect to x are dropped in eq.(l). 
The x dependence of the FS for U c g = is shown in Fig. 
2. We see that the FS becomes round around x w 0.7. 
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FIG. 2. x dependence of to, to, to, as well as the correspond- 
ing FS's. x is expected to approach from x = [/3'-type] to 
x = 1 [/3-type] as the applied pressure increases. 

Here, we calculate the self-energy for the dimer model 
for U e g = 0.4 <~ 0.5eV by using the FLEX approxima- 
tion [5]. Hereafter, we write U e s as U for simplicity. To 
obtain the magnetic transition temperature Tn, we calcu- 
late the Stoner factor without vertex corrections (VC's), 
a s , given by 



as - max { £7 • x°(k, w = 0) }, 



(2) 



where \° (q> ) is given by 

X °(q, u>i) = G(q + k, w, + e„)G(k, e n ), (3) 

k.n 

where G(q + k, u>i + e„) is the Green function given by 
the FLEX approximation, and ui (e n ) is the Matsub- 
ara frequency for boson (fermion). Tn is determined by 
the Stoner criterion; as = 1. In the FLEX approxima- 
tion, however, as does not exceed 1 at finite temper- 
atures in two dimensional systems, which is consistent 
with the Mermin and Wagner theorem. So we deter- 
mine Tn by the condition as = %, where we set as 
(1 - a N ) _1 - 0(100). The AF state will occur through 
the weak coupling between layers, J±. [13] 

Next, We solve the linearized Eliashbcrg equation to 
obtain the SC transition temperature T c . For a singlet- 
pairing case [</>(— k, e n ) — +</>(k, e n )\, it is given by [12] 
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A • 0(k, e„) = VQl - q, e„ - e m ) 



xG(q,e m )G(-q, -e m ) • 0(q,e m ), (4) 

x°(k,^) 



where V(k,w,) = §17 rnT^k^ - 2^ i+c/ x "(k,^) 
T c is given by the condition A = 1 . In the FLEX approx- 
imation, a finite T c is obtained even in two dimensional 
systems irrespective of the Hohenberg theorem. How- 
ever, this approximation gives reasonable T c 's for k-(ET) 
organic compounds and high-T c cuprates [2-4,12]. Note 
that we could not find parameters where the triplet pair- 
ing is dominant in the present model. 

The obtained phase diagram is given in Fig. 3. 64x64 
/c-points and 512 Matsubara frequencies are used. The 
maximum T c under the condition of T c > Tn is about 
16K (18K) at x = 0.8 {x = 0.7) for U = 0.5eV 
(U = 0.4cV). T c decreases as x increases, and T c ~ 6K 
at x = 1 for U = 0.5eV, which is consistent with the ob- 
served T c in several /3-(ET) 2 X compounds. The present 
study also explains the low T c in /3-(ET) 2 X at ambient 
pressure, which corresponds to x = 1 in Fig. 3. On 
the other hand, we could not find the SC phase in the 
metallic region for a; < 0. We comment that T c ~ 4K 
is obtained at x = for U = 0.4eV as shown in Fig. 
3. However, it should be covered with the AF phase be- 
cause of Tn 3> T c . Actually, the Stoner factor at 4K is 
extremely close to 1; as ~ 0.999. 

On the other hand, the Necl temperature for U = 
0.5eV at x = is 28K (20K) under the condition of 
on = 0.99 («n = 0.995), which is also consistent with 
T N = 22K in /?'-(ET) 2 ICl 2 at ambient pressure. The ob- 
tained magnetic order is commensurate [Q = (tt, tt)] for 
x < 0.7 at U = 0.5eV, and it becomes incommensurate 
for x > 0.8 at lower temperatures. The obtained Tn in- 
creases as x departs from 0, which is interpreted as a nat- 
ural consequence of the dimensional increase. It is highly 
demanded to experiment on the pressure dependence of 
Tn- We note that the crossover temperature Tmit deter- 
mined by the condition dp/dT = will be higher than 
Tn because Tmit reflects the Mott transition, which is 
beyond the scope of the FLEX approximation. 

Figure 4 shows the solution of eq.(4) at (k x , k y ) on the 
Fermi surface, where #k = tan~ 1 (k x /k y ). Thus, the ob- 
tained SC order parameter is d a .2_ y 2-wave like, as in the 
high-T c cuprates and the k-(ET) compounds. Note that 
the present x-y coordinate is different from that of the 
original crystal; see Fig. 1. 

Finally, we study the resistivity p and the Hall coef- 
ficients i?H- Considering that the present model lacks 
the four-fold rotational symmetry, we define p and Rn 

as p = ^2/((jI x + Oy~) and i? H = {<Jxy/B z ) ■ p 2 respec- 
tively: They are independent of the choice of the x-y 
coordinate. We calculate them by including the VC's for 
the current to maintain the conservation laws [14]. The 
role of the VC for the current, which is totally dropped 
in the relaxation time approximation, is very important 
in strongly correlated systems. For example, the anoma- 



X°(k,^) 



u. 



lous behaviors of i?u in high-T c cuprates and in k-(ET) 2 X 
compounds are reproduced only when the VC for the cur- 
rent is taken into account adequately [15-17]. 
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FIG. 3. Obtained phase diagram by the FLEX approxima- 
tion. 
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FIG. 4. Obtained FS as well as 0(6» k , iirT) on the FS 
(U = 0.5eV, x = 0.8, T = 15K). 

Figure 5 shows the obtained results for U = 0.5eV. Ex- 
perimentally, p under 8.2GPa continues to increase up to 
300K approximately proportional to T, which is similar 
to p in high-T c cuprates [18]. Such a "bad metal" behav- 
ior of p is well reproduced in the present study. In Fig. 5, 
i?H for x — 0.8 increases below ~ 50K, which is caused by 
the VC's for the current due to the AF fluctuations. How- 
ever, the obtained enhancement of i?H is rather smaller 
than that for high-T c cuprates or /3-(ET) 2 X compounds. 
As for x = 0.9 and 1.0, the temperature dependence of 
i?H is much moderate. It is highly demanded to exper- 
iment on i? H in /3'-(ET) 2 ICl 2 . We stress that the ob- 
served resistivity under high pressure is considered to be 
intrinsic free from the lattice contraction effect, which 
is considerably large in usual organic metals at ambient 
pressure [18]. 
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In more detail, however, it might be better to analyze 
the original tight-binding model in Fig.l, instead of the 
dimer model, because a transport coefficient is in general 
sensitive to the shape of the FS. 
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FIG. 5. Obtained temperature dependences of p and _Rh 
with including the VC's for the current. 

Here, we discuss the validity of the pressure induced 
dimensional crossover in /3'-(ET) 2 ICl 2 . According to 
Taniguchi, the resistivity along the less-conductive direc- 
tion (i.e., c-axis direction) decreases very rapidly as the 
pressure increases, which suggests the system becomes 
2D-likc [18]. In fact, the Q1D nature in /^'-compound is 
ascribed to the accidental cancellation of t\ in the dimer 
model although |i(<?2)| and \t(c)\ are not so small; see Fig. 
1 and Table I. It suggests that the dimensional crossover 
easily occurs. Moreover, /3'-(ET) 2 AuCl 2 remains semi- 
conducting even at 9GPa [18], which means that the /?'- 
type structure in /3'-(ET) 2 AuCl 2 is more robust against 
the pressure. It is natural because Au-Cl bond length is 
about 12% smaller than the I-Cl bond length. 

As a result, the obtained phase diagram, Fig. 3, as 
well as the concept of the SC driven by the dimensional 
crossover, will make sense and be reasonable, although 
the variable x in eq.(l) cannot be interpreted as the pres- 
sure simply. The large density of states (DOS) at the 
Fermi level for x « 0.7 due to the van Hove singularity 
around (±tt,0) (see Fig. 2) accounts for the high T c in 
/?'-(ET) 2 ICl 2 . Actually, we could not obtain higher T c 
( > 5K) under the condition that T c > Tn by the slight 
modification of the original parameters for /3'-(ET) 2 ICl 2 , 
as far as the FS is Q1D. We note that T c in a Q1D sys- 
tem, TMTSF, is very low (w IK), which is recognized by 



the FLEX approximation [19]. 

In the present analysis, we might have underestimated 
the enhancement of the ratio Wb an d/U (ITband being 
the band width) due to the applied pressure. Actually, 
the hydrostatic pressure not only distorts the structure 
(/?' — » f3), but also shortens the cell volume, which en- 
hances each overlap integral. Then, the enhancement of 
t(pl), which is proportional to U e g, is weaker than other 
hopping parameters [11]. 

In summary, we studied the origin of the supercon- 
ductivity in /3'-(ET) 2 ICl 2 based on the FLEX approx- 
imation, by assuming that the structure approaches to 
/3-type under the hydrostatic pressure. Our theory pre- 
dicts that the d-wave superconductivity occurs in /?'- 
(ET) 2 IC1 2 as a result of the 1D-2D dimensional crossover 
owing to the crystal structure change under high pres- 
sure. Moreover, both p and Ru were studied in terms of 
the conserving approximation. Experimental studies on 
the structure of /3'-(ET) 2 ICl 2 as well as the band calcula- 
tions under the condition of high pressure are highly de- 
manded. In addition, the superconductivity in /3-(ET) 2 X 
compound is also understood in the same framework. 

The author is grateful to H. Taniguchi and M. 
Miyashita for valuable discussions on experiments. He is 
also grateful to R. Kato for useful informations on band 
calculation. 
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